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O ■ The problem of choosing the optimal multipath components to be employed at a minimum mean square error 

^__l , (MMSE) selective Rake receiver is considered for an impulse radio ultra-wideband system. First, the optimal finger 

^\j ' selection problem is formulated as an integer programming problem with a non-convex objective function. Then, 
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the objective function is approximated by a convex function and the integer programming problem is solved by 

means of constraint relaxation techniques. The proposed algorithms are suboptimal due to the approximate objective 
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function and the constraint relaxation steps. However, they perform better than the conventional finger selection 
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O ■ algorithm, which is suboptimal since it ignores the correlation between multipath components, and they can get 
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quite close to the optimal scheme that cannot be implemented in practice due to its complexity. In addition to 



the convex relaxation techniques, a genetic algorithm (GA) based approach is proposed, which does not need any 
approximations or integer relaxations. This iterative algorithm is based on the direct evaluation of the objective 
function, and can achieve near-optimal performance with a reasonable number of iterations. Simulation results are 
presented to compare the performance of the proposed finger selection algorithms with that of the conventional 
and the optimal schemes. 
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I. Introduction 

Since the US Federal Communications Commission (FCC) approved the limited use of ultra-wideband 
(UWB) technology [1], communications systems that employ UWB signals have drawn considerable 
attention. A UWB signal is defined to be one that possesses an absolute bandwidth larger than 500MHz or 
a relative bandwidth larger than 20% and can coexist with incumbent systems in the same frequency range 
due to its large spreading factor and low power spectral density. UWB technology holds great promise for 
a variety of applications such as short-range high-speed data transmission and precise location estimation. 

Commonly, impulse radio (IR) systems, which transmit very short pulses with a low duty cycle, are 
employed to implement UWB systems ([2]-[6]). In an IR system, a train of pulses is sent and information 
is usually conveyed by the positions or the amplitudes of the pulses, which correspond to pulse position 
modulation (PPM) and pulse amplitude modulation (PAM), respectively. In order to prevent catastrophic 
collisions among different users and thus provide robustness against multiple-access interference, each 
information symbol is represented by a sequence of pulses, and the positions of the pulses within that 
sequence are determined by a pseudo-random time-hopping (TH) sequence specific to each user [2]. The 
number Nf of pulses representing one information symbol can also be interpreted as pulse combining 
gain. 

Typically, users in an IR-UWB system employ Rake receivers to collect energy from different multipath 
components. A Rake receiver combining all the paths of the incoming signal is called an all-Rake (ARake) 
receiver. Since a UWB signal has a very wide bandwidth, the number of resolvable multipath components 
is usually very large. Hence, an ARake receiver is not implemented in practice due to its complexity. 
However, it serves as a benchmark for the performance of more practical Rake receivers. A feasible 
implementation of multipath diversity combining can be obtained by a selective-Rake (SRake) receiver, 
which combines the M best, out of L, multipath components [7]. Those M best components are determined 
by a finger selection algorithm. For a maximal ratio combining (MRC) Rake receiver, the paths with highest 
signal-to-noise ratios (SNRs) are selected, which is an optimal scheme in the absence of interfering users 
and inter-symbol interference (ISI) [8]-[10]. For a minimum mean square error (MMSE) Rake receiver, 
the "conventional" finger selection algorithm can be defined as choosing the paths with highest signal- 
to-interference-plus-noise ratios (SINRs). This conventional scheme is not necessarily optimal since it 



ignores the correlation of the noise terms at different multipath components. The finger selection problem 
is also studied in the context of CDMA downlink equalization, and recursive sequential search (RSS) 
and heuristic arguments of interference cancellation are proposed to determine finger locations [11]. The 
RSS algorithm selects the fingers one by one in a sequential manner, which reduces the computationally 
complexity significantly; however, it may be quite suboptimal depending on the correlation structure of the 
noise components. The heuristic arguments are employed to determine finger locations when the number 
of fingers is larger than the number of multipath components, which is not applicable in our case due 
to a large number of multipath components in a typical UWB channel. In [12], asymptotic optimality of 
"regular" finger assignments is investigated. However, the behavior of the algorithms for a small number 
of fingers, and performance of other schemes without "regular" assignments are not specified. Finally, the 
finger selection problem is considered in [13] for UWB systems, and a matching pursuit-based technique 
with a quadratic constraint is proposed. However, the performance of the algorithm depends heavily on a 
parameter of the quadratic constraint, which needs to be determined empirically. 

In this paper, unlike the previous approaches, we provide a complete optimization theoretical framework 
for the finger selection problem for MMSE SRake receivers. First, we formulate the optimal MMSE SRake 
as a nonconvex, integer-constrained optimization, in which the aim is to choose the finger locations of 
the receiver so as to maximize the overall SINR. While computing the optimal finger selection is NP- 
hard, we present several relaxation methods to turn the (approximate) problem into convex optimization 
problems that can be very efficiently solved by interior-point methods, which are polynomial-time in the 
worst case, and are very fast in practice. These optimal finger selection relaxations produce significantly 
higher average SINR than the conventional one that ignores the correlations, and represent a numerically 
efficient way to strike a balance between SINR optimality and computational tractability. Moreover, we 
propose a genetic algorithm (GA) based scheme, which performs finger selection by iteratively evaluating 
the overall SINR expression. Using this technique, near-optimal solutions can be obtained in many cases 
with a degree of complexity that is much lower than that of optimal search. 

The remainder of this paper is organized as follows. Section II describes the transmitted and received 
signal models in a multiuser frequency- selective environment. The finger selection problem is formulated 
and the optimal algorithm is described in Section III, followed by a brief description of the conventional 



algorithm in Section IV. In Section V, two convex relaxations of the optimal finger selection algorithm, 
based on an approximate SINR expression and integer constraint relaxation techniques, are proposed. 
The GA-based approach is presented in Section VI as an alternative to the suboptimal schemes based on 
convex relaxations. Simulation results are presented in Section VII, and concluding remarks are made in 
the last section. 

II. Signal Model 

We consider a synchronous, binary phase shift keyed TH-IR system with K users, in which the 
transmitted signal from user k is represented by: 

s -0 = M E d f «&]**(* - jTf - <fT^ (1) 

V * j=—oo 

where pu.{t) is the transmitted UWB pulse, Ef. is the bit energy of user k, Tf is the "frame" time, 
Nf is the number of pulses representing one information symbol, and b[ ', N , G {+1, —1} is the binary 
information symbol transmitted by user k. In order to allow the channel to be shared by many users and 
avoid catastrophic collisions, a TH sequence {<r- }, where c, G {0, 1, ...,N C — 1}, is assigned to each 
user. This TH sequence provides an additional time shift of c: T c seconds to the jth pulse of the fcth 
user where T c is the chip interval and is chosen to satisfy T c < Tf/N c in order to prevent the pulses from 
overlapping. We assume Tf = N C T C without loss of generality. The random polarity codes <r- ' are binary 
random variables taking values ±1 with equal probability [14]- [16]. 

An example TH-IR signal is shown in Figure [T] where six pulses are transmitted for each information 
symbol (Nf = 6) with the TH sequence {2,1,2,3,1,0}. 

Consider the discrete presentation of the channel, ct^ = [a[ ■ ■ ■ a L ] for user k, where L is assumed 
to be the number of multipath components for each user, and T c is the multipath resolution. Note that 
this channel model is quite general in that it can model any channel of the form Xl/=i &i $(t ~ u ) if 
the channel is bandlimited to 1/T C [17]. Then, the received signal can be expressed as 



r ^ = J2\W E T.^^bwN^Ut-jTf-cf^-il-im + aMt), (2) 
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where p rx (t) is the received unit-energy UWB pulse, which is usually modelled as the derivative of p tx (t) 



due to the effects of the receive antenna, and n(t) is zero mean white Gaussian noise with unit spectral 
density. 

We assume that the TH sequence is constrained to the set {0, 1, ... , Nt — 1}, where Nt < N c — L, so 
that there is no inter-frame interference (IFI). However, the proposed algorithms are valid for scenarios 
with IFI as well, and this assumption is made merely to simplify the expressions throughout the paper. 
From the analysis in [18], the results of this paper can easily be extended to the IFI case as well. 

Due to the high resolution of UWB signals, chip-rate and frame-rate sampling are not very practical for 
such systems. In order to have a lower sampling rate, the received signal can be correlated with symbol- 
length template signals that enable symbol-rate sampling of the correlator output [19]. The template signal 
for the Zfh path of the incoming signal can be expressed as 

(i+l)JV)-l 

*£*#) = E dfp^t- 3 T f -cfT c -{l-l)T c l (3) 
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for the ith information symbol, where we consider user 1 as the desired user, without loss of generality. 
In other words, by using a correlator for each multipath component that we want to combine, we can 
have symbol-rate sampling at each branch, as shown in Figure El 

Note that the use of such template signals results in equal gain combining (EGC) of different frame 
components. This may not be optimal under some conditions (see [18] for (sub)optimal schemes). However, 
it is very practical since it facilitates symbol-rate sampling. Since we consider a system that employs 
template signals of the form ©, i.e. EGC of frame components, it is sufficient to consider the problem 
of selection of the optimal paths for just one frame. Hence, we assume Nf = 1 without loss of generality. 

Let C = {li, . . . ,1m} denote the set of multipath components that the receiver collects (Figure |3). 
At each branch, the signal is effectively passed through a matched filter (MF) matched to the related 
template signal in © and sampled once for each symbol. Then, the discrete signal for the /th path can 
be expressed, for the zth information symbol, as 2 

n = sfAbi + n h (4) 

for I = h,..., l M , where A = diagjv^&T, • • • , VEk}, b; = [by ■ ■ ■ b\ K '] T and n L ~ A/"(0 , a 2 .), s^ is 

2 Note that the dependence of n on the index of the information symbol, i, is not shown explicitly. 



a K x 1 vector, which can be expressed as a sum of the desired signal part (SP) and multiple-access 
interference (MAI) terms: 

s, = sr + si MAI) , (5) 

where the kih. elements can be expressed as 



„(SP) 



a, (1) , k = 1 



J k 



and (6) 



k = 2,...,K 



(MAI) 



(i).(fc)v^L ,J*)r( fc ) 



(7) 
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with Jj ^ being the indicator function that is equal to 1 if the mth path of user fc collides with the Zth 
path of user 1, and otherwise. 



III. Problem Formulation and Optimal Solution 

The problem is to choose the optimal set of multipath components, C = {l\, . . . , Im}, that maximizes 
the overall SINR of the system. In other words, we need to choose the best samples from the L received 
samples ri, I = 1, . . . , L, as shown in ©. 

To reformulate this combinatorial problem, we first define an M x L selection matrix X as follows: M 
of the columns of X are the unit vectors ei, . . . , &m (e^ having a 1 at its 2th position and zero elements 
for all other entries), and the other columns are all zero vectors. The column indices of the unit vectors 

determine the subset of the multipath components that are selected. For example, for L = 4 and M = 2, 
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X 



chooses the second and third multipath components. 
10 
Using the selection matrix X, we can express the vector of received samples from M multipath 



components as 

r = XSAb, + Xn, (8) 

where n is the vector of thermal noise components n = [m ■ • ■ni J ] T , and S is the signature matrix given 



by S = [si ■ • • s L ] T , with s t as in ©. 
From ©-Q, © can be expressed as 



r = bf ] v^XaW + XS^ MAI ) Ab, 4 + Xn, (9) 

where S( MAI ) is the MAI part of the signature matrix S. 
Then, the linear MMSE receiver can be expressed as 

k = sign{0 r r}, (10) 

where the MMSE weight vector is given by [20] 

= R- 1 Xc* (1) , (11) 

with R being the correlation matrix of the noise term: 

R = XS( MAI >A 2 (S( MAI )) T X T + a 2 n l. (12) 

The overall SINR of the system can be expressed as 

SINR(X) = %aW) T X T (l + i-XS( MAI )A 2 (S( MAI )) T X T ^| ' Xa' 1 '. (13) 

< \ < J 

Hence, the optimal finger selection problem can be formulated as finding X that solves the following 
problem: 

maximize SINR(X), (14) 

subject to the constraint that X has the previously defined structure. Note that the objective function to be 
maximized is not concave and the optimization variable X takes binary values, with the previously defined 
structure. In other words, two major difficulties arise in solving (fl4l globally: nonconvex optimization 
and integer constraints. Either makes the problem NP-hard. Therefore, it is an intractable optimization 
problem in this general form. 



IV. Conventional Algorithm 

Instead of the solving the problem in (fill) , the "conventional" finger selection algorithm chooses the 
M paths with largest individual SINRs, where the SINR for the /th path can be expressed as 



^ (1) ) 2 

= ( 8 CMAD )r ^(MAD +{J ,. < 15 > 



Ei{a\ 



for / = 1, . . . , L. 



This algorithm is not optimal because it ignores the correlation of the noise components of different 
paths. Therefore, it does not always maximize the overall SINR of the system given in (fT3l) . For example, 
the contribution of two highly correlated strong paths to the overall SINR might be worse than the 
contribution of one strong and one relatively weaker, but uncorrelated, paths. The correlation between the 
multipath components is the result of the MAI from the interfering users in the system. 

V. Relaxations of Optimal Finger Selection 

Since the optimal solution in (fl4ll is quite difficult, we first consider an approximation of the objective 
function in (fT3t . When the eigenvalues of 4-XS( MAI )A 2 (S( MAI )) T X T are considerably smaller than 1, 
which occurs when the MAI is not very strong compared to the thermal noise, we can approximate the 
SINR expression in ( TH51) as follows 3 : 

SINR(X) « ^(c^f X T (l - i-XS( MAI )A 2 (S( MAI )f X^) Xa« (16) 

n \ n / 

which can be expressed as 

SINR(X) « ^-{ (a«) T X T Xa« - ^a«X T XS( MAI )A 2 (S( MAI )) T X T X«« }. (17) 

Note that the approximate SINR expression depends on X only through X T X. Defining x = [x± ■ ■ ■ xl] t 
as the diagonal elements of X T X, x = diag{X T X}, we have xt = 1 if the ith path is selected, and x\ = 

3 More accurate approximations can be obtained by using higher order series expansions for the matrix inverse in <13i . However, the 
solution of the optimization problem does not lend itself to low complexity solutions in those cases. 



otherwise; and ^ i=1 x i = M. Then, we obtain, after some manipulation, 



SINR(x) = ^ jq T x - 4 xTpx j > ( lg ) 

where q = [(af } ) 2 ■ • • (c^T and P = diag{af } • ■ • a^jS^^A^S^^fdiagja^ ■ • ■ af}. 
Then, we can formulate the finger selection problem as follows: 

minimize — x T Px — x T q 

subject to x T l = M, 

Xi6{0,l}, i = l,...,L. (19) 

Note that the objective function is convex since P is positive definite, and that the first constraint 
is linear. However, the integer constraint increases the complexity of the problem. The common way 
to approximate the solution of an integer constraint problem is to use constraint relaxation. Then, the 
optimizer will be a continuous value instead of being binary and the problem (fT9l will be convex. Over the 
past decade, both powerful theory and efficient numerical algorithms have been developed for nonlinear 
convex optimization. It is now recognized that the watershed between "easy" and "difficult" optimization 
problems is not linearity but convexity. For example, the interior-point algorithms for nonlinear convex 
optimization are very fast, both in theory and in practice: they are provably polynomial time in theory, and 
in practice usually take about 25-40 times the computational load of solving a least-squares problem of 
the same dimension to reach a global optimum of a constrained convex optimization problem [21], [22]. 
The interior-point methods solve convex optimization problems with inequality constraints by applying 
Newton's method to a sequence of equality constrained problems, where the Newton's method is a kind 
of descent algorithm with the descent direction given by the Newton step [21]. 

We consider two different relaxation techniques in the following subsections. 
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A. Case-1: Relaxation to Sphere 

Consider the relaxation of the integer constraint in (IT9l) to a sphere that passes through all possible 
integer values. Then, the relaxed problem becomes 

minimize — x T Px — x T q 

subject to x T l = M, 

(2x- l) T (2x- 1) < L. (20) 

Note that the problem becomes a convex quadratically constrained quadratic programming (QCQP) 
problem [21]. Hence it can be solved for global optimality using interior-point algorithms in polynomial 
time. 

B. Case-2: Relaxation to Hypercube 

As an alternative approach, we can relax the integer constraint in (IT9t to a hypercube constraint and 
obtain 

minimize — x T Px — x T q 

subject to x T l = M, 

XG[0,1] L , (21) 

where the hypercube constraint can be expressed as x y and x ■< 1, with y y z meaning that 
2/i > Zi,...,yL > Zl. Note that the problem is now a linearly constrained quadratic programming 
(LCQP) problem, and can be solved by interior-point algorithms [21] for the optimizer x*. 

C. Dual Methods 



We can also consider the dual problems. For the relaxation to the sphere considered in Section IV-A1 
the Lagrangian for d2"0l can be obtained as 

£(x, A, v) = x T ( — P + z/I J x - x T (q - Al + ul) - MX, (22) 
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where X E TZ and v £ 1Z + . 

After some manipulation, the Lagrange dual function, which is the Lagrangian maximized over the 
primal variable x, can be expressed as 

g(X, v) = -i [q + {v - A)1] T f -Lp + vl\ [q + {v - A)l] - MX, (23) 

Then, the dual problem becomes 

minimize -fq + (u - X)1] T I -\p + vl ) [q+ (v - A)l] + MA (24) 

4 ' W / 

subject to v > 0, (25) 

which can be solved for optimal A and v by interior point methods. Or, more simply, the unconstrained 
problem (l24l can be solved using a gradient descent algorithm, and then the optimizer v is mapped to 
v* = max{0, v\. 

After solving for optimal A and /i, the optimizer x* is obtained as 

x* = i (-Lp + i/l) [q+K-A*)l]. (26) 

Note that the dual problem (|2"4"1) has two variables, A and z/, to optimize, compared to L variables, the 
components of x, in the primal problem (l20b . However, an L x L matrix must be inverted for each iteration 
of the optimization of (|2"4"1) . Therefore, the primal problem may be preferred over the dual problem in 
certain cases. 

Similarly, the dual problem for the relaxation in Section IV-BI can be obtained from d2TT) as 

2 

minimize — (q+ fj, - v - Al) T P^(q + n - v - Al) + MX + v T \ (27) 

subject to n,i^h 0. (28) 

It is observed from (l27t that there are 2L + 1 variables and also L x L matrix inversion operations 
for the solution of the dual problem. Therefore, the simpler primal problem (I2TT) is considered in the 
simulations. 
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D. Selection of Finger Locations 

After solving the approximate problem <fT9l> by means of integer relaxation techniques mentioned above, 
the finger location estimates are obtained by the indices of the M largest elements of the optimizer x* . 

Both the approximation of the SINR expression by (fToT) and the integer relaxation steps result in 
the suboptimality of the solution. Therefore, it may not be very close to the optimal solution in some 
cases. However, it is expected to perform better than the conventional algorithm most of the time, since 
it considers the correlation between the multipath components. However, it is not guaranteed that the 
algorithms based on the convex relaxations of optimal finger selection always beat the conventional one. 
Since the conventional algorithm is very easy to implement, we can consider a hybrid algorithm in which 
the final estimate of the convex relaxation algorithm is compared with that of the conventional one, and 
the one that minimizes the exact SINR expression in (IT~3l> is chosen as the final estimate. In this way, the 
resulting hybrid suboptimal algorithm can get closer to the optimal solution. 

VI. Finger Selection Using Genetic Algorithms 

The algorithms in the previous section convert the optimal finger selection problem into a convex 
problem by approximating the SINR expression in dT3l by a concave function in (ll"8l) . and by employing 
the relaxation techniques on the integer constraints. As another way to solve the finger selection problem, 
we propose a GA based approach, which directly uses the exact SINR expression in (fl3l) . and does not 
employ any relaxation techniques. 

A. Genetic Algorithm 

The GA is an iterative technique for searching for the global optimum of a cost function [23]. The 
name comes from the fact that the algorithm models the natural selection and survival of the fittest [24]. 

The GA starts with a population of chromosomes, where each chromosome is represented by a binary 
string 5 . Let N ipop denote the number of chromosomes in this population. Then, the fittest N pop of these 
chromosomes are selected, according to a fitness function. After that, the fittest iVg 00 d chromosomes, 
which are also called the "parents", are selected and paired among themselves {pairing step). From each 

'Although we consider only the binary GA, continuous parameter GAs are also available [23]. 
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chromosome pair, two new chromosomes are generated, which is called the mating step. In other words, 
the new population consists of N good parent chromosomes and N good children generated from the parents 
by mating. After the mating step, the mutation stage follows, where some chromosomes (the fittest one 
in the population can be excluded) are chosen randomly and are slightly modified; that is, some bits in 
the selected binary string are flipped. After that, the pairing, mating and mutation steps are repeated until 
a threshold criterion is met. 

The GA has been applied to a variety of problems in different areas [23]-[25]. Also, it has recently been 
employed in the multiuser detection problem [26]-[28]. The main characteristics of the GA algorithm is 
that it can get close to the optimal solution with low complexity, if the steps of the algorithm are designed 
appropriately. 

B. Finger Selection via the GA 

In order to be able to employ the GA for the finger selection problem we need to consider how to 
represent the chromosomes, and how to implement the steps of the iterative optimization scheme. 

A natural way to represent a chromosome is to consider the "assignment vector" x in (fT8t . which 
denotes the assignments of the multipath components to the M fingers of the RAKE receiver. In other 
words, %i = 1 if the ifh path is selected, and Xj = otherwise; and J2i=i x i = M. 

Also, the fitness function that should be maximized can be the SINR expression given by (fT3t . Note 
that, given a value of x, SINR(X) can be uniquely evaluated. By choosing this fitness function, the fittest 
chromosomes of the population correspond to the assignment vectors with the largest SINR values. 

Now the pairing, mating and mutation steps need to be designed for the finger selection problem: 

1 ) Pairing: The assignments to be paired among themselves are chosen according to a weighted random 
pairing scheme [23], where each assignment is chosen with a probability that is proportional to its SINR 
value. In this way, the assignments with large SINR values have a greater chance of being chosen as the 
parents for the new assignments. 

2) Mating: From each assignment pair, the two new pairs are generated in the following manner: Let 
xx and x 2 denote two finger assignments, and let p Xl and p X2 consist of the indices of the multipath 
components chosen as the RAKE fingers. Then, the indices of the new assignments are chosen randomly 
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from the vector p = [p Xl p X2 ]. If the new assignment is the same as x x or x 2 , the procedure is repeated 
for that assignment. For example, consider a case where L — 10 and M = 4. If x x = [1 1 1 10 0] 
(p Xi = [1 4 7 8]) and x 2 = [0 1 1 1 1 0] (p X2 = [2469]), then the new assignments are chosen 
randomly from the set p = [14 7 8 2 4 6 9]. For example, the new assignments (children) could be 
x 3 = [1 1 1 1 0] (p X3 = [1 2 4 9]) and x 4 = [0 1 1 1 1 0] (p X4 = [467 9]). 

Note that by designing such a mating algorithm, we make sure that a multipath component that is 
selected by both parents has a larger probability of being selected by the new assignment than a multipath 
component that is selected by only one parent does. 

3) Mutation: In the mutation step, an assignment, except the best one (the one with the highest SINR), 
is randomly selected, and one 1 and one of that assignment are randomly chosen and flipped. This 
mutation operation can be repeated a number of times for each iteration. The number of mutations can 
be determined beforehand, or it might be defined as a random variable. 

Now, we can summarize our GA-based finger selection scheme as follows: 

• Generate N ipop different assignments randomly. 

• Select iVp p of them with the largest SINR values. 

. Pairing: Pair N good of the finger assignments according to the weighted random scheme. 

• Mating: Generate two new assignments from each pair. 

• Mutation: Change the finger locations of some assignments randomly except for the best assignment. 

• Choose the assignment with the highest SINR if the threshold criterion is met; go to the pairing step 
otherwise. 

In the simulations, we stop the algorithm after a certain number of iterations. In other words, the 
threshold criterion is that the number of iterations exceeds a given value. As the number of iterations 
increases, the performance of the algorithm increases, as well. The other parameters that determine the 
tradeoff between complexity and performance are A^ ipop , A'pop, N good , and the number of mutations at 
each iteration. 
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VII. Simulation Results 

Simulations have been performed to evaluate the performance of various finger selection algorithms 
for an IR-UWB system with N c = 20 and Nf — 1. In these simulations, there are five equal energy 
users in the system (K = 5) and the users' TH and polarity codes are randomly generated. We model the 
channel coefficients as a L = sign(o!;)|Q;;| for I = 1, ...,L, where sign(a ? ) is ±1 with equal probability 
and \ai\ is distributed lognormally as £A/"(/i;, er 2 ). Also the energy of the taps is exponentially decaying 
as E{|ai| 2 } = fi e" A( ' _1) , where A is the decay factor and Ef=i E {M 2 } = l ( s ° ft = (1 - e ~ A )/( 1 ~ 
e~ AL )). For the channel parameters, we choose A = 0.1, a 2 = 0.5 and hi can be calculated from /i; = 
0.5 ln( 1 1 J e _ AI , ) — X(l — 1) — 2<t 2 , for / = 1, . . . , L. We average the overall SINR of the system over 
different realizations of channel coefficients, TH and polarity codes of the users. 

In Figure |3J we plot the average SINR of the system for different noise variances when M = 5 
fingers are to be chosen out of L = 15 multipath components. For the GA, A^p = 32, A pop = 16, and 
-^good = 8 are used, and 8 mutations are performed at each iteration. As is observed from the figure, 
the convex relaxations of optimal finger selection and the GA based scheme perform considerably better 
than the conventional scheme, and the GA get very close to the optimal exhaustive search scheme after 
10 iterations. Note that the gain achieved by using the proposed algorithms over the conventional one 
increases as the thermal noise decreases. This is because when the thermal noise becomes less significant, 
the MAI becomes dominant, and the conventional technique gets worse since it ignores the correlation 
between the MAI noise terms when choosing the fingers. 

Next, we plot the SINR of the proposed suboptimal and conventional techniques for different numbers 
of fingers in Figure 0] where there are 50 multipath components and Et/N = 20. The number of chips 
per frame, N c , is set to 75, and all other parameters are kept the same as before. In this case, the optimal 
algorithm takes a very long time to simulate since it needs to perform exhaustive search over many different 
finger combinations and therefore it was not implemented. The improvement using convex relaxations of 
optimal finger selection over the conventional technique decreases as M increases since the channel is 
exponentially decaying and most of the significant multipath components are already combined by all 
the algorithms. Also, the GA based scheme performs very close to the suboptimal schemes using convex 
relaxations after 10 iterations with iVi pop = 128, A pop = 64, A goo d = 32, and 32 mutations. 
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Finally, we consider an MAI-limited scenario, in which there are 10 users with E\ — 1 and Ef. — 10 
\/k ^ 1, and all the parameters are as in the previous case. Then, as shown in Figure |5J the improvement by 
using the suboptimal finger selection algorithms increase significantly. The main reason for this is that the 
suboptimal algorithms consider (approximately) the correlation caused by MAI whereas the conventional 
scheme simply ignores it. 

VIII. Concluding Remarks 

Optimal and suboptimal finger selection algorithms for MMSE-SRake receivers in an IR-UWB system 
have been considered. Since UWB systems have large numbers of multipath components, only a subset 
of those components can be used due to complexity constraints. Therefore, the selection of the optimal 
subset of multipath components is important for the performance of the receiver. We have shown that the 
optimal solution to this finger selection problem requires exhaustive search which becomes prohibitive 
for UWB systems. Therefore, we have proposed approximate solutions of the problem based on Taylor 
series approximation and integer constraint relaxations. Using two different integer relaxation approaches, 
we have introduced two convex relaxations of the optimal finger selection algorithm. 

Moreover, we have proposed a GA based iterative finger selection scheme, which depends on the direct 
evaluation of the objective function. In each iteration, the set of possible finger assignments is updated in 
search of the best assignment according to the proposed GA stages. 

The three contributions of this paper are the formulation of the optimization problem, the convex 
relaxations, and the GA based scheme. In the first, the formulation is globally optimal, but the solution 
methods for this non-convex nonlinear integer constrained optimization must use heuristics due to its 
prohibitive computational complexity. In the second, the suboptimal schemes focus on relaxations where 
globally optimal solutions are guaranteed, but the problem statement is relaxed. In the third, the problem 
statement remains the same as the original, intractable NP hard one, but the solution method is by local 
search heuristics. 
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Fig. 1. An example time-hopping impulse radio signal with pulse-based polarity randomization, where Nf — 6, N c = 4, the time hopping 
sequence is {2,1,2,3,1,0} and the polarity codes are {+1,+ 1,-1,+1,-1,+ 1}. 
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Fig. 2. The receiver structure. There are M multipath components, which are combined by the MMSE combiner. 



20 



30 



25 



Optimal 

Conventional 

Sphere 

Hypercube 

Genetic Algorithm, 10 Iteration 



CD 
DC 



CD 
CD 

CO 

i— 

CD 
> 
< 



20 




E b /N (dB) 



Fig. 3. Average SINR versus E^/Nq for M = 5 fingers, where Ej, is the bit energy. The channel has L = 15 multipath 
components and the taps are exponentially decaying. The IR-UWB system has N c = 20 chips per frame and Nf = 1 frame 
per symbol. There are 5 equal energy users in the system and random TH and polarity codes are used. 
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Fig. 4. Average SINR versus number of fingers M, for Eb/No = 20dB, N c = 75 and L = 50. All the other parameters are 
the same as those for Figure [3] 



22 




Number of Fingers 



Fig. 5. Average SINR versus number of fingers M. There are 10 users with each interferer having lOdB more power than 
the desired user. All the other parameters are the same as those for Figure |4] 



